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Observations show that the Cosmic Microwave Background (CMB) contains tiny variations at 
the 10 level around its black-body equilibrium temperature. The detection of these temperature 
fluctuations provides to modern Cosmology evidence for the existence of primordial density per- 
turbations that seeded all the structures presently observed. The vast majority of the cosmological 
information is contained in the 2-point temperature function, which measures the angular corre- 
lation of these temperature fluctuations distributed on the celestial sphere. Here we study such 
angular correlations using a recently introduced statistic-geometrical method. Moreover, we use 
Monte Carlo simulated CMB temperature maps to show the equivalence of this method with the 
2-point temperature function (best known as the 2-Point Angular Correlation Function). We also in- 
vestigate here the robustness of this new method under possible divisions of the original catalog-data 
in sub-catalogs. Finally, we show some applications of this new method to simple cases. 



INTRODUCTION 

According to modern Cosmology, in the primordial 
Universe ionized barionic matter and photons were 
tightly coupled until the recombination epoch when the 
hydrogen and helium nuclei hold their electrons to form 
neutral atoms, leaving photons to freely stream toward 
us. These photons, nowadays observed as a faint resid- 
ual relic radiation termed the Cosmic Microwave Back- 
ground Radiation (CMB), have a well described black- 
body Planckian spectrum as the FIRAS instrument on 
the COBE satellite showed in 1992 (see e.g. and ref- 
erences therein). 

Besides this thermal equilibrium feature, the CMB pre- 
served evidences of small variations in temperature (or 
intensity) from one part of the microwave sky to an- 
other, for this called CMB temperature fluctuations (also 
known as CMB anisotropics in the literature), which are 
interpreted as being originated during their interaction 
with the primordial matter-density fluctuations. Since 
at that time the density fluctuations were correlated at 
specific angular scales characteristic of the physical pro- 
cesses involved (e.g., 0, 0), such correlations were im- 
printed in the angular distribution of the the CMB tem- 
perature fluctuations. This fact turns of fundamental im- 
portance the study of the CMB angular correlations, in 
order to understand the evolution of the matter-density 
fluctuations to form the structures -like galaxies and 
galaxy clusters- observed today. 

In the last decades, especially after the COBE mis- 
sion, observational cosmology evolved astonishingly fast, 
in such a way that there has been an increasing improve- 
ment in the quality of measuring the CMB temperature 
fluctuations. While with COBE data the CMB studies 
were restricted to large angular-scales, with the highly 
precise and excellent angular resolution of WMAP data, 
analysis at all angular-scales are now possible. 



In what follows we present the basics of a recently pro- 
posed method suitable for the study of the angular cor- 
relations in the sky distribution of cosmic objects. Our 
main interest is focussed in CMB data like those released 
by the WMAP satellite Q, for this we shall apply the 
method to CMB temperature fluctuations maps (briefly 
CMB maps) , and also investigate the robustness of the 
method under possible partitions of the original data in 
sub-catalogs. Finally, we shall show the equivalence of 
our method with the 2-Point Angular Correlation Func- 
tion (2-PACF), and this will be done both in a theoretical 
way using the definitions of these approaches as well as in 
a computational way using Monte Carlo simulated CMB 
maps. 

THE PASH METHOD 

Recently, Bernui & Villela 0, introduced a new 
method to scrutiny the angular correlations in the dis- 
tribution of cosmic objects in the celestial sphere. 
Given a catalog with n cosmic objects (sharing common 
physical properties), the method consists in first calcu- 
late the n(n— 1)/2 angular distances of all pairs-of-objects 
and then construct the normalized histogram of the num- 
ber of pairs with a given angular separation versus their 
angular separation. The Expected-PASH (EPASH) is 
the histogram obtained for the ideal distribution case. 
The difference between the PASH, caculated from the 
observational data, and the EPASH can reveal signifi- 
cant departures of the former histogram with respect to 
the later, if the signal-to-noise ratio is sufficiently large. 
This geometric-statistical tool is called the Pair Angular 
Separation Histogram (PASH) method. 

Catalogs with a large number of objects, as high- 
resolution CMB maps, can be suitable divided in a set 
of comparable sub-catalogs, provided that they contain 
a similar number of objects sharing analogous physical 
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properties. After the partition of the original catalog 
in, say K, sub-catalogs one performs a PASH with each 
one of these sub-catalogs and average them to obtain 
the Mean-PASH (MPASH). As we shall see below, the 
MPASH function is independent of the number of sub- 
maps K in which the original map is divided, provided 
that the minimum number is K m i n = 2: one sub-map for 
the positive CMB and the other for the negative ones. 
Finally we perform, and plot, the difference between the 
MPASH and the EPASH. 

The method has the advantage that it does not de- 
pend on cosmological parameters or data other than the 
angular position in the celestial sphere of the cosmic ob- 
jects listed in the catalog to be analised. The method 
is also applicable to incomplete sky maps, including dis- 
connected regions, on condition that the EPASH was ob- 
tained under similar characteristics (i.e. the same num- 
ber of objects per catalog, and where objects are located 
in the same patch of the celestial sphere). Notice that 
the WMAP CMB maps do not include the monopole and 
dipole contributions, therefore for such data the EPASH 
has to obtained taking into account these features. 

We firstly describes the basics of the method by show- 
ing how to obtain a PASH given a catalog of data. A 
catalog C is a list of n cosmic objects containing suitable 
information about their angular coordinates on the celes- 
tial sphere and some physical properties; here we assume 
that the objects of a given catalog have common physical 
properties. 

We divide the interval (0, 180°] in Nuns bins of equal 
length Sj = 180° /Nbi ns , where each sub-interval has 
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with center in r y i = (i — i) 5j. Next, we calculate the 
N = n(n — l)/2 angular distances corresponding to the 
distances between all pairs-of-objects. We denote by (p(j) 
the number of pairs of objects in C separated by a dis- 
tance 7 G (0, 180°]. Then, 
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is the normalized counting of the number of pairs of ob- 
jects, i.e. X^i'" ^7 = 1? separated by an angle n f i , 
that lies in the sub-interval J^. 

Let p, q any two objects p, q G C, with angular co- 
ordinates {9 p ,<p p ) 1 {6 q ,4>q), respectively, where 9 p , 9 q G 



and 



G [0,360°]. The unit vectors 



[0,180°j cm 1U Vp , 
(i.e., \rip\ = \n q \ = 1) denote the position of the two ob- 
jects p, q on the celestial sphere. We denote by p{n p ) the 
probability densities of the objects distributed on the ce- 
lestial sphere <S 2 and listed in C . Then, the probability 
density that two objects p, q G C be separated an angular 
distance 7 G [0, 180°] is defined by 

^(7)= // dn p dQ q p(n p )p(n q )6(j-d(ri p ,ri q )), (2) 



where S is the Dirac-delta, and d(n p , n q 
arccos [cos 9 p cos 8 q + sin 9 p sin 9 q cos(4> p — (j> q )]. 
Hence, the PASH is simply defined by 
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However, if the interval Ji is small enough (as shall be 
considered here) , a suitable approximation for the PASH 
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Notice that p(n p )dfl p is the probability of the object p to 
be in the sky patch dfl p G S 2 . Of course, the probability 
of finding p in the whole S 2 should be equal to 1, 



p{n p ) dfl p 



1 . 



Thus, a purely isotropic and normalized density dis- 
tribution of objects observed in the full-sky S 2 , reads 
p{n p ) — j^. With this information, a direct calculation 
in eq.J2J) gives, 
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is the EPASH for the case of isotropically distributed 
objects on the full-sky S 2 . 

Notice that a PASH V^) satisfies the normalisation 

property: J 180 V^d-y — 1, which let us to perform 
the mean of an arbitrary number of histograms, for sub- 
catalogs containing comparable number of objects. This 
is a very useful condition since produces a sort of nor- 
malisation of our method in the sense that area of the 
function MPASH-minus-EPASH is zero: 
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THE 2-POINT ANGULAR CORRELATION 
FUNCTION 

The Two-Point Angular Correlation Function (2-PACF) 
is defined by 

C( 7 )^(T(n p )T(n 9 )}, (7) 

where Tip j Tig cirC such that cos 7 = n p ■ n qi and 
T(ri p ), T(n q ) are the temperature fluctuations of the pix- 
els p, q of a given CMB map, respectively. Notice that, in 
some literature (see, e.g. 0) this definition appears with 
T(ri p )/To instead of T(n p ), where To is the black-body 
equilibrium temperature of the CMB. 
Since rii , i — p, q is a vector with angular coordinates 
(9i,(j)i), i — p,q, one can expand T(n^) in the spherical 
harmonics, so 



T(fii) = y^^ai m Yi 
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where ag m are the multipole moments, with zero mean 
(o-im) = for Gaussian temperature fluctuations. In a 
statistically isotropic Universe the variance of the multi- 
pole moments is independent of m, which means that we 
can define 
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and the set of Ci is termed the angular power spectrum 
of the CMB map. 

Using eqs. (JSJ and (J0J, and well-known properties of 
the spherical harmonics in the definition @ , the 2-PACF 
for a statistically isotropic Universe can be written as 



C(l) = ^^2(2£+l)C e P,(co S1 ), 



(10) 



where Pi are the Legendre polynomial of oreder I. 
define the function 



We 
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and call it the normalized 2-PACF because the area 
between the curve and the horizontal axis is zero, as 
can be verified from the bottom-plot appearing in Fig- 
ure 1. Moreover, assuming a Harrison- Zeldovich (HZ) 
scale-invariant spectrum and for large angular scales (i.e., 
I < 100) we have Cg = Q/(^(£ +1)1 where Q is the 
quadrupole normalisation constant []], |a] . In Figure 1 we 
plotted both functions, C( 7 ) and C"°™( 7 ), for the case 
of a statistically isotropic Universe with HS spectrum. 

When we use the formula 1)10(1 , instead of the definition 
eq. 0, a crucial point appears in the method of the 
2-PACF. In fact, intrinsically to the method there is a 
large uncertainty in the evaluation of the Ces, called the 
Cosmic Variance, and is given by SCg/Ci = \j2j {21 + 1). 
This reveals no more than our ignorance in determining 
the real values of the low-order Ce (i.e. for £ = 2, 3, 4, ...). 



THE 2-PACF VERSUS THE PASH-minus-EPASH 
FUNCTION 

Here we shall prove the equivalence between the PASH 
method and the 2-PACF. 

The above definition of the 2-PACF, eq. (JJJ, can also 
be written as (see Q, pag. 195, eq. (5.28), where £(x) = 
£(cos 7 )^C( 7 )) 

C( 7 ) + 1 

= A dtt q p(n q ) dSlp p{n p )S (cos 6 P - cos( 7 -I- q ))5{<f> p -<f> q ), 



where A is a constant, and S is the Dirac-delta. Multi- 
plying both sides by which according to eq. [3] is the 
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FIG. 1: Here we consider a statistically isotropic Universe 
with HZ spectrum. Moreover, motived in WMAP data where 
the monopole and dipole contributions are not included, here 
we consider that the sum index starts with I = 2, i.e. with the 
quadrupole moment. Top: We plotted the 2-PACF, given by 

C (7) = %Y.i=2 W+Tj p c( cosr y)- Bottom: We plotted the 
normalized 2-PACF, which defined by C norm ( 7 ) = ^ C( 7 ); 
this 2-PACF is called normalized because the area between 
the curve and the horizontal axis is zero. 



EPASH for the full-sky catalog, and using properties of 
the Dirac-delta we obtain 
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which means that when 



then 7 = 9 P 



other words the only contribution to the integral comes 
from the equality 7 = d{n pi n q ). This implies that the 
right-hand-side of this equation is equal to the definition 
of the PASH, eq. J3J), except for the constant A. Thus, 
choosing suitably the value of A, we formally obtain the 
equivalence between both approaches 
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FIG. 2: The red-points (including the standard devia- 
tion bars at la) correspond to average of 1000 MPASH- 
minus-EPASH functions from Monte Carlo CMB maps gen- 
erated under the statistical isotropy assumption, and the 
black-line corresponds to the expected 2-PACF C norm (j) = 
^ fff+i"; fi( cos 7) f° r a statistically isotropic Uni- 

verse. We realize that the excellent coincidence between both 
results is a validation of the equivalence between the MPASH- 
minus-EPASH function and the normalized 2-PACF. 

Moreover, we also performed here a computational ver- 
ification of this equivalence. For this, we carry out the av- 
erage of 1000 MPASH-minus-EPASH functions from an 
equal number of Monte Carlo CMB full-sky maps gen- 
erated under the statistical isotropy and HZ spectrum 
assumptions. Our results, observed in Figure 2, show 
the excellent coincidence between both functions which 
is a computational verification of the equivalence between 
the MPASH-minus-EPASH function and the normalized 
2-PACF. 

Therefore, we conclude that the two approaches, 
the MPASH-minus-EPASH and the normalized 2-PACF 
(C norm (7) = ^Cfr)) are equivalent. However, it is 
worth to mention here that the normalized version of the 
2-PACF seems to be more suitable than the original 2- 
PACF version, because it leds to scrutiny the angular 
correlations at all angular scales. As observed in the top- 
plot of Figure 1, the original 2-PACF is not well-defined 
for small angular scales [0, ~ 20°]. 

ROBUSTNESS OF THE MPASH CALCULATION 

Here we shall study the robustness in the obtention of 
the MPASH, in the case that it could be obtained con- 
sidering different partitions (or ways to generate the sub- 
catalogs) of the original catalog (e.g. the temperature 
data in a CMB map). We shall deal with this problem 
by directly showing our results of calculating the MPASH 



in different partitions (different number of sub-catalogs 
K) of a given WMAP CMB temperature map. 

A CMB full-sky map is the celestial sphere partitioned 
in a set of equal- area pixels, where to each pixel is as- 
signed a temperature value; so the above mentioned phys- 
ical property is the CMB temperature fluctuation value 
of each pixel. To obtain the angular correlations of a 
CMB map, we order the set of pixels according to increas- 
ing values of the pixels temperature. Next, we divide this 
ordered set of pixels in K disjoint sub-sets (termed sub- 
catalogs) , all of them with a comparable number of pix- 
els. Afterward, we compute the MPASH averaging the 
K PASHs calculated from each of the K sub-catalogs. 
As we shall show in Figure 6, the MPASH function is 
independent of the number of sub-catalogs K in which 
the original map is divided, provided that the minimum 
number is K m { n = 2: one sub-map for the positive CMB 
and the other for the negative ones. However, it is clear 
that we can divide the original data in several ways, al- 
ways following these criteria. 

Here we consider the CMB map obtained using WMAP 
data, specifically the co-added Q+V+W map tj, which 
is a weighted combination of the 8 high frequency 
differencing assemblies (DA): Ql, Q2, VI, V2, Wl, W2, 
W3, and W4 (see: http:\\lambda.gsfc.nasa.gov\product 
\map\current\IMaps _cleaned.cfm), that is obtained ac- 
cording to 

_ Ej=3 r -r/°o,J 
1 co-added v-^io -, / 2 ' 

where Tj is the CMB map for the DA /, and Cg j is the 
noise per observation for the DA /. 

To investigate whether the MPASH function is in- 
dependent of the number of sub-catalogs K in which 
the original catalog is divided due to a large amount 
of pixels, we study first the angular correlations in the 
North-Galactic and South-Galactic spherical cap regions 
of 70° of aberture of the co-added WMAP map. Our 
neat results shown in Figure 3, for the analysis of seven 
different partitions of the above mentioned regions of 
the co-added CMB map, demonstrates that the MPASH- 
minus-EPASH function (and of course also the MPASH) 
is actually independent of K. 

Finally, in Figure 4 we show the same result, but 
this time for the full-sky WMAP CMB map, termed the 
'cleaned' map which was produced by combining all 
the individual frequency WMAP maps using a Wiener fil- 
ter algorithm in order to reduce or eliminate foregounds 
(specially those coming from our Galaxy). 

SIMPLE APPLICATIONS OF THE PASH 
METHOD 

In order to illustrate the MPASH-minus-EPASH func- 
tion revealing signatures in simple cases, here we con- 
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FIG. 3: MPASH-minus-EPASH functions for CMB data in 
the 70° North and South Galactic spherical caps of the co- 
added WMAP map The seven curves plotted in each 
figure (top and bottom) correspond to the cases where the 
number of histograms was: K = 4, 16, 17, 20, 25, 40, and 80, 
respectively. The K = 4 case corresponds to the degraded 
map N a ide = 64, and all the others to N 3 ide = 128 As 
observed, the coincidence of this rainbow of curves proves that 
the MPASH-minus-EPASH is independent of the number K. 

sider three CMB maps where the angular correlations 
are known a priori: 

• a statistically isotropic CMB map, that is there are 
no angular correlations between the temperature 
fluctuations of the map; the results of our analysis 
are shown in Figure 5. 

• a pure dipole map; the results of our analysis are 
shown in Figure 6. 

• a pure quadrupole map; the results of our analysis 
are shown in Figure 7. 
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FIG. 4: MPASH-minus-EPASH functions for the full-sky 
WMAP 'cleaned' map [To| . The four curves plotted here 
correspond to the cases where the number of histograms to 
construct the MPASH was: K = 40, 81, 100, and 130, respec- 
tively. As observed, all the coloured-curves coincide showing 
again that the MPASH-minus-EPASH function is indepen- 
dent of the number K. 
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FIG. 5: Statistically isotropic CMB map and the MPASH- 
minus-EPASH analysis for these data. The MPASH-minus- 
EPASH function is just a statistical (noisy) oscillations 
around the horizontal axis, revealing just the absence of an- 
gular coorelations in the data. 
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FIG. 6: A CMB Dipole map and the MPASH-minus- EPASH 
analysis for these data. The MPASH-minus-EPASH function 
results in a Sine function, trivially revealing a tipical dipole- 
signature (with some statistical noise). 
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